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Abstract 

The process of radiation from high-energy electron in oriented single crys¬ 
tal is considered using the method which permits inseparable consideration 
of both coherent and incoherent mechanisms of photon emission. The total 
intensity of radiation is calculated. The theory, where the energy loss of 
projectile has to be taken into account, agrees quite satisfactory with avail¬ 
able CERN data. It is shown that the influence of multiple scattering on 
radiation process is suppressed due to action of crystal field. 
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Recently authors developed a new approach to analysis of pair creation by 
a photon in oriented crystals This approach not only permits to consider 
simultaneously both the coherent and incoherent mechanisms of pair creation by 
a photon but also gives insight on the Landau-Pomeranchuk-Migdal (LPM) effect 
(influence of multiple scattering) on the considered mechanism of pair creation. In 
the approach the polarization tensor of photon was used which includes influence of 
both external held and multiple scattering of electrons and positrons in a medium 
j2]. In the present paper the analysis of process of radiation from a high-energy 
electron in oriented crystal includes inhuence of both an external held and the 
multiple scattering of electron. This makes possible indivisible consideration of 
both coherent and incoherent mechanisms of photon emission as well as analysis of 
inhuence of the LPM ehect on radiation process. 

The properties of radiation are connected directly with details of motion of 
emitting particle. The momentum transfer from a particle to a crystal we present 
in a form q =< q > -|-qs, where < q > is the mean value of momentum transfer 
calculated with averaging over thermal(zero) vibrations of atoms in a crystal. The 
motion of particle in an averaged potential of crystal, which corresponds to the 
momentum transfer < q >, determines the coherent mechanism of radiation. The 
term q^ is attributed to the random collisions of particle which dehne the incoherent 
radiation. Such random collisions we will call ’’scattering” since < q* >= 0. If the 
radiation formation length is large with respect to distances between atoms forming 
the axis, the additional averaging over the atom position should be performed. 

Under some generic assumptions the general theory of the coherent radiation 
mechanism was developed in [3]. If the electron angle of incidence 'do (the angle 
between electron momentum p and the axis (or plane)) is small 'do -C Vo/m, where 
Vo is the characteristic scale of the potential, the held E of the axis (or plane) 
can be considered constant over the pair formation length and the constant-held 
approximation (magnetic bremsstrahlung limit) is valid. In this case the behavior 
of radiation probability is determined by the parameter 

£ E 

X = 

m Eq 

where £ is the electron energy, m is the electron mass, Eq = m? je = 1.32-10^® V/cm 
is the critical held, the system h = c = 1 is used. The very important feature of 
coherent radiation mechanism is the strong enhancement of its probability at high 
energies (from factor ~ 10 for main axes in crystals of heavy elements like tungsten 
to factor ~ 170 for diamond) comparing with the Bethe-Heitler mechanism which 
takes place in an amorphous medium. If 'do 3> Vo/m the theory passes over to 
the coherent bremsstrahlung theory (see 13,0 0) Side by side with coherent 
mechanism the incoherent mechanism of radiation is acting. In oriented crystal this 
mechanism changes also with respect to an amorphous medium j7]. The details 
of theory and description of experimental study of radiation which confirms the 
mentioned enhancement can be found in 0. The study of radiation in oriented 
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crystals is continuing and new experiments are performed recently isi, i. 

At high energies the multiple scattering of radiating electron (the LPM effect) 
suppresses radiation probability when e > Se- In an amorphous medium (or in 
crystal in the case of random orientation) the characteristic electron energy starting 
from which the LPM effect becomes essential is Se ~ 2.5 TeV for heavy elements USI 
and this value is inversely proportional to the density. In the vicinity of crystalline 
axis (just this region gives the crucial contribution to the Bethe-Heitler mechanism) 
the local density of atoms is much higher than average one and for heavy elements 
and at low temperature the gain could attain factor ~ 10^. So in this situation the 
characteristic electron energy can be Eq ~ 2.5 GeV and this energy is significantly 
larger than ’’threshold” energy St starting from which the probability of coherent 
radiation exceeds the incoherent one. It should be noted that the main contribution 
into the multiple scattering gives the small distance from axis where the held of 
crystalline axis attains the maximal value. For the same reason the LPM effect in 
oriented crystals originates in the presence of crystal held and nonseparable from 
it. This means that in problem under consideration we have both the dense matter 
with strong multiple scattering and high held of crystalline axis. 

Below we consider case VQ/m. Than the distance of an electron from axis 
Q as well as the transverse held of the axis can be considered as constant over the 
formation length. For an axial orientation of crystal the ratio of the atom density 
n(p) in the vicinity of an axis to the mean atom density na is 


where 


n{x) 

na 




Vi 


^0 


e(o)’ 


Xo = 


TldUnal ’ 


2ul 


Q 

at 


( 2 ) 

(3) 


Here g is the distance from axis, ui is the amplitude of thermal vibration, d is the 
mean distance between atoms forming the axis, is the ehective screening radius 
of the axis potential (see Eq.(9.13) in jH]) 


U{x) 




In 1 + 


X + T] ^ 


In 1 + 


Xo + V, 


(4) 


The local value of parameter x{x) (see Eq.(^) which determines the radiation 
probability in the held Eq.(jD) is 


X{x) 


dU{g) e 2y^ VqE _ e 

dg rrfi [x-\-r]){x + t] + 1)' rrfias Eg 


(5) 


The parameters of the axial potential for the ordinarily used crystals are given 
in Table 9.1 in p. The particular calculation below will be done for tungsten 
crystals studied in [S]. The relevant parameters are given in Table 1. It is useful 
to compare the characteristic energy Eq with ’’threshold” energy Et for which the 
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radiation intensity in the axis field becomes equal to the Bethe-Maximon one. Since 
the maximal value of parameter x(a;): 


Xm = = t(\/l + 16i;(l + rfj - 1 - 2))), = — 


( 6 ) 


is small for such electron energy (st -C Em), one can use the decomposition of 
radiation intensity over powers of y (see Eq.(4.52) in jH]) and carry out averaging 
over X. Retaining three terms of decomposition we get 


jF ^ (ao{7^) - ai{7])xs + a2{v)xl + •••), 


ao{r]) = (1 + 2r]) In - 2 , 




ai{ri) = 
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TT 


64 


1 


02 (h) = 64 


(1 + 277) 


Vv V^ + v 
1 


4 (\/l + h - 


M^ + v) 


+ 30 - 12(1+ 577 ( 1 + 77 )) In 


1+77 


(7) 


The intensity of incoherent radiation in low energy region e < Et £m (see 
Eq.(21.16) in [H] and Eq. (jA.18jl in Appendix A) 


am^ E 

I = - - 9o 


47r Ee 


9o = ^ + -j- 


1 + 34.4 (x^ In X + 2 . 54 x 2 ^ 


uT 


- h{^ 

18 


> / = Jf{x)e ’'1—, 


( 8 ) 


where 


TTt 1 

= mZVXin^L„ ' + 2 - 

1 111/“^/^ 

h{z) = --[l + {l + z)e^Ei{-z)], a= ^ , 


/(O = Re [7/^(1 + 7O - V'(l)] = 




n{n^ + ^2) ’ 


(9) 


here 'ip{z) is the logarithmic derivative of the gamma function, Ei(z) is the integral 
exponential function, /(^) is the Coulomb correction. For x = 0 fhis intensity 
differs from the Bethe-Maximon intensity only by the term h{u\/a?‘) which reflects 
the nongomogeneity of atom distribution in crystal. For tti -C a one has h{u\/a^) ~ 
— (1 + C)/2 + ln(a/ 77 i), C = 0.577.. and so this term characterizes the new value of 
upper boundary of impact parameters Ui contributing to the value < > instead 

of screening radius a in an amorphous medium. 
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Conserving in Eq. m only the main (the hrst) term of decomposition, which 
corresponds to the classical radiation intensity, neglecting the corrections in Eq. (jH|) 
(^fQ = 1, X = 0), using the estimate Vq ~ Zajd and Eqs.Q, Q, we get 


2Tiao{r]) 


63 


Lod 

aoiv) 


MeV, 
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where the distance d is taken in units 10“® cm. Values of et found using this 
estimate for tungsten, axis < 111 >, d=2.74 -10“® cm are consistent with points of 
intersection of coherent and incoherent intensities in Fig.l (see Table 1). For some 
usable crystals (axis < 111 >, room temperature) one has from Eq. |Tn|) 


ei(C(d)) ~ 0.47 GeV, ^^(Si) ~ 2.0 GeV, et{Ge) ~ 1.7 GeV, (11) 


so this values of et are somewhat larger than in tungsten except the diamond very 
specihc crystal where value of et is close to tungsten one. 

For large values of the parameter Xm S> £m) the incoherent radiation intensity 
is suppressed due to the action of the axis held. In this case the local intensity of 
radiation can by written as (see Eq.(7.129) in P) 


jinc ^ 29r(l/3) £ 
3V62430 £e 

Here we take into account that 
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1 (0727+ti^ 
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Averaging the function (x(x))“^/^ and \nx{x){x{x))~‘^^^ over x according with 
Eq.(jHl) one can End the effective value of upper boundary of the transverse mo¬ 
mentum transfer (oc rnxH^ instead of m) which contributes to the value < >. 

Using the obtained results we determine the effective logarithm L by means of 
interpolation procedure 

L = Log, ^ = ^0 + ^ In (l + 70xV) ■ (14) 

Let us introduce the local characteristic energy (see Eq. O) 

edx) = (15) 

^{x)9 9 

In this notations the contribution of multiple scattering into the local intensity for 
small values of Xm and e/eo has a form (see Eq.(15) in [TT] I 
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Integrating this expression over x with the weight I/xq we get 


jLPM 
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( 17 ) 


It should be noted that found Eq. (ini) has a good accuracy only for energy much 
smaller (at least on one order of magnitude) than Eq (see discussion after Eq.(15) 

in dH). 

The spectral probability of radiation under the simultaneous action of multiple 
scattering and an external constant field was derived in [H] (see Eqs.(7.89) and 
(7.90)). Multiplying the expression by a; and integrating over uj one obtains the 
total intensity of radiation I. For further analysis and numerical calculation it is 
convenient to carry out some transformations 


1. Changing of variables: u —>• az//2, r —>• 2t/a, (z/r — vt). 

2. Turn the contour of integration over t at the angle —7r/4. 
One finds after substitution t —>■ \/2t 


„ i xo oo 

, am r ydy f ax ^ , f ^ , tt 

I{e) = / -- / —Gr{x,y), Gr{x,y) = / Fr{x,y,t)dt - 

/ttjI — j/jxq J 4 


Fr{x,y,t) + ihr)(p2{t) F ^ , h 


4x^(x) 


^ 9 9 ’ 


y = -, M = -^, ipi{t) = {i-l)t + br{l + i){(p2{t)-t), 

£ 1-y 

= -^tanh-^, =- '~~ 1 = -5 

z/o V2 sinh(y2z/of) 


(18) 


where 

r2 = l + {l-yf, r3 = 2(l-|/), (19) 

y £c[x) 

oj is the photon energy, the function £c{x) is defined in Eq. (US and x(x) is defined 
in Eq. ©■ The expression for the spectral probability of radiation used in the above 
derivation can be found from the spectral form of Eq.(16) in P {dW/dy = ujdW/de) 
using the standard QED substitution rules: £ —»• — uj —>• —a;, e'^de —*• uj‘^duj and 
exchange uJc{x) —>• 4£c(x). 

The inverse radiation length in tungsten crystal (axis < 111 >) 1/L‘^'’(e:) = 
I{e)/£ Eq. (fTK|l . well as coherent contribution 1/L^{£) = I^{£)/£ Eq.(j2SI) and in¬ 
coherent contribution 1/L*"''^(e:) = F'^^[e')le Eg. ^Tj) are shown in Fig.l for two 
temperatures T=100 K and T=293 K as a function of incident electron energy e. 
In low energy region (e < 0.3 GeV) the asymptotic expressions Eqs.((7j) and (jHI) are 
valid. One can see that at temperature T=293 K the intensity 1^[e) is equal to 
at £ ~ 0.4 GeV and temperature T=100 K the intensity I^{£) is equal to 
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P"'^{e) at e ~ 0.7 GeV. The same estimates follow from comparison of Eqs. dZD and 
Q; see also Eq. m- At higher energies the intensity dominates while the 

intensity decreases monotonically. 

The inverse radiation length given in Fig.l can be compared with data directly 
only if the crystal thickness I <C L^'^{e) (thin target). Otherwise one has to take 
into account the energy loss. The corresponding analysis is simplihed essentially 
if I < = (max(/(£)/£))“^. The radiation length L'^'^{e) varies slowly on the 

electron trajectory for such thicknesses. This is because of weak dependence of 
L^'^{e) on energy in the region L^'^{e) ~ imin relatively large value of 

L^^{e) ^ i^rnin region where this dependence is essential but variation of 

energy on the thickness I is small. For W, axis < 111 >, T=293 K one has 
j^min _ 220 at energy £ = 300 GeV, see Fig.l. For this situation dispersion can 
be neglected (see discussion in Sec. 17.5 of jH]) and energy loss equation acquires 
the form 


1 de 


= -L^ 




e dl ' ' e 

In the hrst approximation the hnal energy of electron is 


( 20 ) 


E, = Eoexp (-llL"{£a)), 


( 21 ) 


where £o is the initial energy. In the next approximation one has 
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( 22 ) 


If the dependence of L^''{e) ^ on £ is enough smooth it’s possible to substitute the 
function by an average value with the weight l/e\ 
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Numerical test conhrms this simplihed procedure. Using it we hnd 
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Enhancement of radiation length (the ratio of Bethe-Maximon radiation length 
LBM ) in tungsten, axis < 111 >, T=293 K is shown in Fig.2. The curve 

1 is for the target with thickness I = 200 /rm, where the energy loss was taken 
into account according using the simplihed procedure Eq. (|2l. The curve 2 is for a 
considerably more thinner target, where one can neglect the energy loss. The only 
available data are from [Hj. The measurement of radiation from more thin targets 
is of evident interest. 

In order to single out the inhuence of the multiple scattering (the LPM ehect) on 
the process under consideration, we should consider both the coherent and incoher¬ 
ent contributions. The probability of coherent radiation is the hrst term (uq = 0) 
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of the decomposition of Eq. (HHD over z/q. The coherent intensity of radiation is 
(compare with Eq.(17.7) in [H]) 


XQ 

I^ie) = I Hx) 
0 


dx 

Xo' 


(25) 


Here I{x) is the radiation intensity in constant held (magnetic bremsstrahlung 
limit, see Eqs. (4.50), (4.51) in jH]). It is convenient to use the following represen¬ 
tation for I{x) 

2 A+200 / 2 \ s 7 

= / (^) r(l-s)r(3;,-l)(2i,-l)(<.^-;, + 2)-^, 

ZTT J \ O I COS TTS 

A—zoo 

^ < A < 1. (26) 

The intensity of incoherent radiation is the second term (oc of the mentioned 
decomposition. In Appendix A the new representation of this intensity is derived, 
which is suitable for both analytical and numerical calculation: 


/'“(.) = (27) 

bUvr Eq j Xq 

where J(x) is dehned in Ea. f)A.16fl . 

The contribution of the LPM effect in the total intensity of radiation I Eq. m 
is dehned as 

jLPM ^ j_jF _ Jinc ^28) 

The relative contribution (negative since the LPM ehect suppresses the radiation 
process) A = is shown in Fig.3. This contribution has the maximum 

A ~ 0.8% at £ ~ 0.7 GeV for T=293 K and A ~ 0.9% at e ~ 0.3 GeV for 

T=100 K or, in general, at e ~ The left part of the curves is described quite 

satisfactory by Eq. m- For explanation of right part of the curves let us remind 
that at e ^ Em the behavior of the radiation intensity at a; ~ is dehned by 
the ratio of the contributions to the momentum transfer of multiple scattering and 
that of the external held on the formation length If (see Eq.(21.3) in [Hj) 


k = 
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< Hs > 

< q >2 


{wlff 


^0 


.y-4/3 ^ 
A,m 


a am? 2/3 - 1/3 

- Xrk = -—X m 5 


/ 




Fr. 


£0 \ / 


(29) 


where w is an acceleration in an external held. The linear over k term determines 
the contribution into intensity of incoherent process: 1/U'^^{e S> Em) k/L^{E) ~ 











am ?/The LPM effect is defined by the next term of decomposition over 
k (oc and decreases with energy even faster than Moreover one has to 

take into account that aX e> Eg the contribution of relevant region x ~ 771 into the 
total radiation intensity is small and 1/L^{e) decreases with the energy growth as 
For such energies the main contribution gives the region x ~ = [e/EsY'^'^ 

and 1/L^'^{e) increases until energy e ~ (see Fig.l). This results in essential 
reduction of relative contribution of the LPM effect A. 

It’s instructive to compare the LPM effect in oriented crystal for radiation 
and pair creation processes. The manifestation of the LPM effect is essentially 
different because of existence of threshold in pair creation process. The threshold 
energy ujm is relatively high (in W, axis < 111 >, ujm ~ 8 GeV for T=100 K and 
ujm ~ 14 GeV for T=293 K). Below 0 ;^ influence of field of axis is weak and the 
relative contribution of the LPM effect attains 5.5 % for T=100 K pp. There is 
no threshold in radiation process and becomes larger than at much lower 
energy Et and starting from this energy the influence of field of axis suppresses 
strongly the LPM effect. So the energy interval in which the LPM effect could 
appear is much narrower than for pair creation and its relative contribution is less 
than 1 % in W, axis < 111 >. Since value of Et depends weakly on Z fEa. (fTn|) i. 
Em OC Z~^ (Eqs.Q. Q) and Eq oc Z~‘^ (Eq.®) the relative contribution of the 
LPM effect A for light elements significantly smaller. Thus, the above analysis 
shows that influence of multiple scattering on basic electromagnetic processes in 
oriented crystal (radiation and pair creation) is very limited especially for radiation 
process. 
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A Appendix 

New representation of the intensity of the incoherent radia¬ 
tion in external field, asymptotic expansions 

In the expression for the intensity of incoherent radiation enters following inte¬ 
gral over photon energy uj (see Eq.(21.21) in P): 


J{X) = J [y^Mz) + Mz)) + 2(1 

0 


y)f 2 {z) dy, z = 


y 


2/3 


M^-y). 


(A.i) 


where y = uj/e, the functions fi{z) and f 2 {z) are defined in the just mentioned 
equation in P: 


hiz) = z^T{z)-3z^r{z)-z^ 

f 2 {z) = (/ 32;)T(^) - 5z‘^T'{z) - z^, (A.2) 


here T(^) is the Hardy function: 


T{z) = J sin ^zt + — ^ dt 


Introducing the variable ?/ = |//(x(l — //)) we obtain 
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JM = / 


:(/. + A) + 


{i+vxy 


(i + hx)‘ 


dy 


{l + yxf 




where 


A,2(x) = J fiA^) 
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Integrating Ea. ()A.5fl by parts we find 

/i,2(oo) 


dy 


(1 + vx) 




JiAx) = 
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AM 


y^/^dy 
(1 + hX)’ 


(A.3) 


(A.4) 


(A.5) 


(A.6) 


Since the integral Ea. dA.bfl for separate terms of functions /(. 2 (^) diverges, one has 
to transform it to an another form. We represent the functions f[ 2 {z) in terms of 
derivative of the Hardy functions 


f[{z)=zM\z)-3zM, fM = M^\z) - 5 zM + 3M, (A.7) 
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where we used equations 

zT{z) = r'{z) + l, + + (A.8) 

Now we will show that 

OO r> 

o 


/ f'i, 2 i.z)d'n = -J z^!’^ f[^^{z)dz = 0. 


(A.9) 


Using Eq. (US) and integration by parts, one can reduce all the integrals in Eq. ()A.9jl 
to the form 


^^,,.dz ^ f ciz r I iT'^ \ , 


OO OO 

dz 


0 0 
OO \ 2 

2 


= ^Re \ I dx\ = 0 . 


(A.IO) 


The last equation permits one to rewrite Ea. dA.hfl as 

JiAx) = + A,2(x), iiAx) =xj 


z^dz 


1 + X^3/2 


(A.ll) 


Entering in Ea. IjA.lUl expression (1 + u) ^ we present as contour integral 

^ . A+ioo 

“ • ./o . . „ 


(1 + u) 2 J sinvrs 

A —200 


ds, u = —1 < A < 0. 


Using the standard form of the Hardy function one has 


d 


= —Im / exp i zr + — dr = Im / (ir)^ exp i rr + — dr 


dz 

0 ' ' ' ' 0 ' ' 

(A.13) 

Substituting in the integral in Ea. IjA.lljl the functions f[ 2 {z) in the form Eqs. (US), 
()A.13jl and integrating over the variables r and r we obtain 


A+200 


= -^ J 


X 


^ 1 , 2 ( 5 ) 


ds 


A —200 


a/3/ r (1 + s/2) sin^(7rs/2) ’ 


(A. 14) 


where r(s) is the gamma function, 


3s 


Ai(s) — r I ——h 6 1 + 3 r I ——h 51, 


3s 


3s 


Ms) = r— + 6+5r— + 5+3r— + 4. 


3s 


3s 


(A.15) 
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Substituting Eas. ()A.14ll . (lA.lbfl into Ea. ()A.ll|l and using Ea. ()A.4ll . we get after 
change of variable s —> 2s, displacement of integration contour and reduction of 
similar terms the hnal expression for J{x) 


A-tt 


X^*r(l + 3s) ds 1 

V Vi ^ --<A<0 (A.16) 

1 (sj sm TTS o 

A —200 


where 


R{s) = 15 + 43s + 31s^ + 283“ + 12s‘' 


(A.17) 


In the case x "C 1, closing the integration contour on the right, one can calculate 
the asymptotic series in powers of y 


J(X) = 15 + 516x2 (In ^ - C 1 + 1893x' + ... ~ 15 


1 - 34.4x^ In-2.542 


X 


(A.18) 


In the case x S> 1 it is convenient to present the integral Ea. ()A.16jl in the form 




f ^r(l-s)r(l + 3s)fl(s)—, -t<A<0 (A.19) 

2 J 3^ sm TTS 3 


A —200 


Closing the integration contour on the left one obtains the series over the inverse 
powers of x 


587rr(l/3) 6287r3^/®r(2/3) 

“ 81 ■ 3V6;^2/3 + 243x^/3 

, 1887rr(l/3) 
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Figure captions 

Fig.l The inverse radiation length in tungsten, axis < 111 > at different tem¬ 
peratures T vs the electron initial energy. Curves 1 and 4 are the total effect: 

= I{e)/e Eq.((THl) for T=293 K and T=100 K correspondingly, the curves 
2 and 5 give the coherent contribution I^[e)/e Eq.(j2SI), the curves 3 and 6 give the 
incoherent contribution Eq. ()27j) at corresponding temperatures T. 

Fig.2 

Enhancement (the ratio in tungsten, axis < 111 >, T=293 K. The 

curve 1 is for the target with thickness I = 200 /im, where the energy loss was 
taken into account (according with Eq. (j24jl ). The curve 2 is for a considerably 
more thinner target, where one can neglect the energy loss (L^-f —> The data 

are from |H]. 

Fig.3 

The relative contribution of the LPM effect A (per cent) in tungsten, axis 
< 111 >. Curve 1 is for T=100 K and curve 2 is for T=293 K. 
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Table 1 Parameters of radiation process of the tungsten crystal, axis < 111 > 

for two temperatures T 


T(K) 

K)(eV) 

Xo 

Vi 

V 

eo(GeV) 

et(GeV) 

^.(GeV) 

£m(GeV) 

h 

293 

413 

39.7 

0.108 

0.115 

7.43 

0.76 

34.8 

14.4 

0.348 

100 

355 

35.7 

0.0401 

0.0313 

3.06 

0.35 

43.1 

8.10 

0.612 
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